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Abstract 



u 
o 

_C We consider upper and lower bounds for maxmin allocations of a 

completely divisible good in both competitive and cooperative strate- 
gic contexts. These bounds are based on the convexity properties of 
the range of utility vectors associated to all possible divisions of the 

(T*) good. We then derive a subgradient algorithm to compute the exact 

value up to any fixed degree of precision. 

^ 1 Introduction 

d 

The notion of what is fair in the allocation of one or more infinitely divisible 
goods to a finite number of agents with their own preferences has long been 
debated. Predictably, no agreement has been reached on the subject. The 
situation is often exemplified with children (players) at a birthday party 
who are around a table waiting for their slice of the cake to be served, with 
the help of some parent (an impartial referee). If we think about a special 
class of resolute children who are able to specify their preferences in terms 
of utility set functions, the parent in charge of the division could ease his 
task by using a social welfare function to summarize the children's utility 
values. Among the many proposals, the maxmin - or Rawlsian - division 
was extensively studied in the seminal work of Dubins and Spanier TO] . 
who showed the existence of maxmin optimal partitions of the cake for any 
completely divisible cake and their main properties. They also showed that 
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when a condition of mutual appreciation holds (assumption (MAC) below) 
any optimal partition is also equitable, i.e., it assigns the same level of utility 
for each child. 

The study of the maxmin optimal partitions and their properties has 
continued in more recent years. In particular, its relationship with other 
important notions such as efficiency (or Pareto optimality) and, above all, 
envy-freeness has been investigated with alternating success: each maxmin 
partition is efficient, but while for the two children case Brams and Taylor 
[B] showed that it is also envy-free, the same may not hold when three or 
more children are to be served, as shown in Dall'Aglio and Hill [9]. 

It is worth pointing out the relationship with n player bargaining solu- 
tions. If we think about the division as deriving from a bargaining procedure 
among children, it is straightforward to show that the bargaining solution 
proposed by Kalai [12], in the case where all the players' utilities are nor- 
malized to 1, coincides with the equitable maxmin division. Therefore, if the 
conditions proposed by Dubins and Spanier hold, the two solutions actually 
coincide. 

Little attention has been devoted, however, to finding optimal maxmin 
partitions with one notable exceptions: the case of two players with additive 
and linear utility over several goods has been considered by Brams and Taylor 
[6J, with the Adjusted Winner procedure. 

For the case of general preferences (expressed as probability measures, i.e. 
nonnegative and countably additive set functions normalized to 1) and arbi- 
trary number of players, Legut and Wilczinski [T7] gave a characterization of 
the optimal maxmin allocation in terms of weighted density functions. More- 
over, Elton et al. pj] and Legut [H] provided lower bounds on the maxmin 
value. The optimization problem was later analysed by Dall'Aglio [7j. The 
general problem was reformulated as the minimization of a convex function 
with a particular attention to the case where the maxmin allocation is not 
equitable and the allocation of the cake occurs in stages to subsets of players. 
No detail, however, was given on how to proceed with the minimization. 

In most of the fair division literature, little is assumed about the strate- 
gic behaviour of the children. Brams and Taylor [6] discuss the issue of the 
manipulability of the preferences: in most cases children may benefit from 
declaring false preferences. A different approach takes into account the pos- 
sibility for the children to form coalitions after (Legut [15] and Legut et al. 
[16J) or before (Dall'Aglio et al. [8]) the division of the cake. In both cases 
coalitional games are defined and analysed. In the case of early cooperation 
among children, the game is based on a maxmin allocation problem among 
coalitions, each one having a joint utility function and a weight. The first 
properties of the game are studied in [8]. It turns out that the analysis of 
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the game is made harder by the difficulties in computing the characteristic 
function, i.e., the value associated to each coalition. The tools we introduce, 
therefore, become essential in computing such values, as well as any synthetic 
value, such as the Shapley value, associated to the game. 

The coalitional maxmin problem is indeed a generalization of the classical 
maxmin problem introduced by Dubins and Spanier. Therefore, we consider 
a common approach to set up an algorithm which, at each step, will compute 
an approximating allocation, together with lower and upper bounds for the 
maxmin value. The algorithm is based on a subgradient method proposed 
by Shor [19] and it yields an approximation of the optimal allocation with 
any fixed degree of precision. 

In Section 2 we describe the maxmin fair division problem with coali- 
tions through the strategic model of interaction among players in [S] and 
the geometrical setting employed in [TJ, [2] and [7]. In Section 3 we present 
the upper and the lower bounds for the objective value. In Section 4 we fit 
the Subgradient Method to our problem and we derive a procedure where 
the optimal value and the optimal partition are computed up to a desired 
precision and we provide a numerical example where we describe two fair di- 
vision games and we compute the corresponding Shapley values. Some final 
considerations are given in Section 5. 

2 The model and the maxmin fair division 
problem with coalitions 

We represent our completely divisible good as the set C, a Borel subset 
of M. n , and we denote as 13(C) the Borel a— algebra of subsets of C. Let 
iV = {1, . . . , n) be the set of players, whose preferences on the good are 
Hi, . . . , \i n , where each fih z G iV, is a probability measures on (C,B(C)). 
By the Radon- Nikodym theorem, if v is a non-negative finite- valued measure 
with respect to which each fii is absolutely continuous (for instance we may 
consider v = ^ 6iV A*i), then, for each A G B(C), 

Hi(A)= [ fcdv VzGiV, 

J A 

where fi is the Radon-Nikodym derivative of fn with respect to v. 
We will consider the following assumptions: 

a) complete divisibility of the good (CD): For each i G N and each A G 13(C) 
such that Hi(A) > 0, there exists a measurable set B C A such that 
Hi(A D B) > and jii(A n B c ) > 0. 
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b) mutual absolutely continuity (MAC): If there exists % £ N and A £ B(C) 

such that fii(A) > 0, then for each j ^ % fij(A) > 0. 

c) relative disagreement (RD): For each pair i,j £ N and each A £ £>(C) 

such that Hi(A) > and fij(A) > 0, there exists a measurable set 
Bci such that ^ + 

Throughout the rest of the work we will assume that (CD) always holds, 
while (MAC) and (RD) are useful, though restrictive, assumptions that we 
will employ only when strictly needed. 

For any h £ N, let (Ai, . . . , Ah) be an /i-partition, i.e., a partition of 
the good C into h measurable sets. Let Tlh be the class of all /i-partitions. 
How do players behave in the division procedure? In the simplest case, each 
player competes with the others to get a part of the cake with no strategic 
interaction with other players. Each (A\, . . . , A n ) £ IT n determines a division 
of the good in which player i, i £ N gets the share Ai with value ^(Ai). Here, 
individual players seek an allocation with values as high as possible. A fair 
compromise between the conflicting interests is given by maxmin allocation 
(A\, . . . , A* n ) £ n n that achieves 



v m := max < min \ . (1) 
04i,... ) A„)en„ ^ hen n J 

Here v m denotes the maxmin value in the classical fair division problem. 
With a completely divisible good, the allocation (A*, . . . , A*) is fair (or 
proportional), i.e. Hi(A*) > - for all i £ N. Moreover, if (MAC) holds, it is 
also egalitarian, i.e. fii(A*) = [ij(A*) for all i,j £ iV (see [ID])- Therefore, 
under this assumption, an optimal allocation is also the bargaining solution 
proposed by Kalai and Smorodinsky [T3] (see also Kalai [T2]). 

Dall'Aglio et al. [S] proposed a strategic model of interaction, where play- 
ers, before the division takes place, gather into mutually disjoint coalitions. 
Within each coalition, players pursue an efficient allocation of their collective 
share of the cake. 

Let Q be the family of all partitions of N and, for each r £ Q, let |T| = m, 
m < n, and let M — {1, . . . , m} be the coalitions indexes set. Thus, players 
cluster into coalitions specified by the partition T = {Si, . . . , S m }. For each 
j £ M and each coalition Sj, players state their joint preferences as follows 



{ D i}ieS i partition of B J B 



£>(A)= / 



VsAB) = max Y]fi i (D i )= / f s Ax)dx (2) 



with fsj(x) = max ieSj . f { (x), B £ B(C) and {Di} ieSj £ B(C). The utility 
fiSj(B) of coalition Sj will be divided among its members in a way that 
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prevents any of them to break the coalition in search of a better deal. Once 
the global coalition structure is known, a fair allocation of the cake among 
the competing coalitions is sought. In this context, assigning the same value 
to all coalitions could yield an unfair outcome. Fairness here must consider 
the different importance that coalitions may assume and this is taken into 
account by a weight function w : 2 N — > M. + . 

In this framework, each coalition takes the role of a single player in Equa- 
tion 0. Following Kalai [12], when coalitions in T are formed and the weight 
function w is considered, players should agree on a division of the cake which 
achieves the following value 



Each coalition can evaluate its performance in the division by considering 
the following coalitional game 



where T$ = {S, {j}j^s}- The value f](S, w) can be interpreted as the minimal 
utility that coalition 5* is going to receive in the division when the system of 
weight w is enforced, independently of the behaviour of the other players. 

A crucial question lies in the definition of the weight system. We consider 
two proposals: 

• w car d = 15*1, S C N. This is certainly the most intuitive setting. Al- 
though very natural, this proposal suffers from a serious drawback, 
since players participating in the game rj(-,w) may be better off wait- 
ing to seek for cooperation well after the cake has been divided (see 



• «V e = fj, s (U ie sA*), S C N, where (A\, A%, . . . , A*) is the partition 
maximizing ([!]). By seeking early agreements among them, players 
will be better off than postponing such agreements until the cake is 
cut. The above mentioned problem is overcome at the cost of a less 
intuitive (and more computationally challenging) formulation (see [8]). 
It is interesting to note that to find these weights we need to solve ([I]). 

It is easy to verify that, for each S C N, 




(3) 



r)(S,w) = w(S)v(T s ,w) 



S C N 



(4) 



0); 



Tj(S, W ca rd) < V(S, W prc ) 

with equality if S = N or S = {i} where % G N. 



5 



The optimization problem ([3]) can be seen as an infinite dimensional as- 
signment problem. In principle we could attribute any point of the cake C to 
any of the participating players (provided certain measurability assumptions 
are met). For very special instances this becomes a linear program: when 
the preferences have piecewise constant densities, or when the cake is made 
of a finite number of completely divisible and homogeneous parts. 

The fully competitive value special instance of the cooperative 

case, since v m = v(Ti,Wi), with I\ = {{1}, . . . , {n}} and wi({i}) = 1 for 
each % G N. Therefore, we focus on the cooperative case alone. 

2.1 A geometrical setting 

We now describe a geometrical setting already employed in [TJ, [2J, [7] and [H] 
to explore fair division problems. In what follows we consider the weighted 
preferences and densities, fij and /j", given respectively by 

w _ t~^Sj /, w _ fSj 

j ~ w{Sj) j ~ w(Sj) ' 

The partition range, also known as Individual Pieces Set (IPS) (see [2]) is 
defined as 

V := {K(£i), • • -,lC(B m )) : (B x , . . . , B m ) G U m } C R™. 

Let us consider some of its features. Each point p G V is the image, un- 
der (/ii, . . . , /i n ), of an m-partition of C. Moreover, V is compact and, if 
(CD) holds, V is also convex (see [18J). Therefore, v(T,w) = max {x > : 
(x, x, . . . , x) n V 7^ 0}. So, the point v(T, w) is the intersection between the 
Pareto frontier of V and the egalitarian line 

£ = {x G R m : x x = x 2 = . . . = x m }. (5) 

3 Upper and lower bounds for the maxmin 
value 

We turn our attention to a simpler optimization problem that may have an 
unfair solution, but it provides easy-to-compute upper and lower bounds for 
the original problem. These bounds depend on a weighted maxsum partition, 
which we can derive through a straightforward extension of a result by Dubins 
and Spanier [10J. Let A m _i denote the unit (m — 1)— simplex. 
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Proposition 3.1. (see ITty Theorem 2], J7| Proposition 4-3]) Let a £ A m „i 
an<i /e£ B a = (-Bf , . . . , B^) be an m— partition of C . If 

akfk(x) > a h fj?(x) for all h,k £ M and for all x £ (6) 

m 

(Sf,...,S«)e argmax ^a^J^). (7) 

[Bi,...,B m )ETl m 

The value of this maxsum problem is itself an upper bound for problem 
([3]). For each choice of a £ A m _i, we have a maxsum partition B a = 
(Bf, . . . , B^) corresponding to a. 

Definition 3.2. The partition value vector (PVV) u a = (mj,...,m^) is 
defined by 

Uj = /jLj(Bj), for each j — 1, . . . , m. 

The PVV u a is a point where the hyperplane X]j=i a j x j = k touches the 
partition range V, so u a lies on the Pareto border of V. Moreover, for any 
a £ A m _i there exists at least one PVV (see [2]). We are ready to state the 
first approximation result. 

Proposition 3.3. Let g : A m _i — > 1R + be as follows: 

g(a) := / max {otjf™(x)}dx. 
Jc J' eM 

Then, 

v(T,w) < q(a) < max/. 

V 1 _ M\ / _ j( _ M j 

Proof. Following [71 Proposition 4.3] we know that the hyperplane that touches 
V at the point / is defined by the equation 

^aiXi = g(a) 

Since a £ A m _i, this hyperplane intersects the egalitarian line £ defined in ^ 
at the point (g(a), . . . ,g(a)). Since the hyperplane is located above V, this 
point lies above the maxmin point with coordinates (v(T,w), . . . ,v(T,w)). 
Therefore 

g(a) > v(T,w) 

Finally, since g(a) is a weighted average of the values («",..., u^), it 
follows that g(a) < max,, gA / w". 

□ 
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The function g was already considered in [7j, where it was shown that g 
is convex, and v(T,w) = min a6 A m _i g(oi)- 

We now turn our attention to a lower bound for v(T,w). Although we 
will see later only one PVV is enough to assure such a bound, we give a 
general result for the case where several PVVs have already been computed. 
We derive the second approximation result through a convex combination of 
these easily computable points in V, which lie close to v(T, w). The following 
result generalizes Theorem 3 in [13] and Theorem 1.1 in [TTj . 



Proposition 3.4. Let u = (u\, . . . ,u m ) a partition value vector such that 



Uh = max Uj. 

j=l,...,m 



Then, 



v(T,w) > v(u) := — - — — > min-u,. (9) 

1 -I- l^jjLh ^(C) J 



Proof. Let us consider the following vectors 

e* = (0,...,0,^(C),0,...,0) qeM q ? h, (10) 

where fJ^(C) is the weighted joint utility of the whole cake by coalition S q . 
Now, consider the convex hull of the PVV u and the m—1 points e q , q ^ h, 

V := {t h u + V 9 : (th, ■■■,t m )e A m _i} 

q+h 

The lower bound we are looking for is the intersection point between V and 
the egalitarian line £ from ^ (see Figure [T]). Let us denote this point as 
(x w , . . . , x w ). Without loss of generality, let us suppose h — 1. Then, we 
obtain (x w , . . . , x w ) as follows: 

tiUi + + . . . + = x w 
hu 2 + t 2 (J%{C) + . . . + = x w 

t x u m + + . . . + t m /j^(C) = x w 
ti + 1 2 + . . . + t m = 1 

We are dealing with a linear system with m + 1 unknown quantities, 

t\ , t 2 , • • • , t mi X w . 

Thus, by Cramer's rule, we get JC' n_; 3jS 
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•1 U r 







/C(C)/ 



— ti^. 



"1 



rig^l (C) 1 + Zlg^l V J%(C) 



1 I (Ul-Uq) 

1 "I" Z^ g ^i ^(c) 



where the second equality derives by suitable exchanges of rows and columns 
in the denominator matrix. In fact, we get the second one after an even 
number of exchanges on the first: m successive exchanges of the last row 
until it reach the first position, and m successive exchanges of the last column 
until it reach the first position. So the two matrices in the denominator have 
the same determinant. It is easy to verify that ti > 0, for every i e N. 

Finally, since the lower bound belongs to the convex hull of the PVV u 
and the m — 1 vectors e q , q ^ h, it is not less than the minimum component 
of each vector, in particular v(u) > mmj^M Uj- 

□ 

An illustration of the position of the bounds with respect to the partition 
range in the case of two coalitions is shown in Figure [TJ 



4 The subgradient method 

In the previous section we have seen that for each choice of the coefficients 
vector a we can derive upper and lower bounds for v(T,w). We describe a 
way of improving the coefficients a so that eventually the bounds will shrink 
to the desired value. 

Since in general g is a non-differentiable convex function, we can rely 
on a simple minimizing algorithm developed by Shor |19j . the subgradient 



9 




Figure 1: Upper and lower bounds for the two-coalition case. 



method. In particular, since the domain of g is constrained, we must consider 
an extension, the projected subgradient method, which solves constrained 
convex optimization problems. Let us start by describing the method through 
some basic definitions and the essential convergence result. 

Definition 4.1. Let D be a closed convex set and let \\ ■ || be the Euclidean 
norm. The projection of z G IR n on D is denoted by p(z) and it is defined as 

p(z) — argmin \\x — z\\. (11) 

x£D 

Definition 4.2. Let f be a convex function with domain D and let xq an 

interior point of D. A vector 7(20) called a subgradient or a generalized 
gradient of f at xq if it satisfies 

f(x)-f(xo)>(j(x ),x-x ) forallxeD. (12) 

Moreover, 7 is a bounded subgradient of f if there exists G G M. + such 
that ||7(x)|| < G for all x G D. 

We denote as d x f(x) the set of subgradients of a convex function / at 
any interior point x of the / domain. 

Definition 4.3. A sequence {s t }^ of positive numbers is called diminishing 
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step size rule if it satisfies the conditions: 



lim s t = 0, (13) 

t— i> + 00 

+00 

J2 s t = +oo. (14) 

t=0 

The subgradient method minimizes a non-differentiable convex function 
which has a bounded set of minimum points and at least one bounded sub- 
gradient. This procedure returns a minimum value for the function moving 
a point in the domain in the opposite direction of a bounded subgradient by 
a step belonging to a diminishing step size rule (see [19]). If the domain of 
the function is constrained, then the point is projected in the domain (see 
[3]). We recall the general result 

Proposition 4.4. (see |3J/, j Let f be a convex function defined on D C 
M. m , which has a bounded set of minimum points D* and let j(x) G d x f(x) 
be a bounded subgradient. Moreover, let {s t }fj^ be a diminishing step size 
rule. Then for any x° G D the sequence {x 1 }'^ generated according to the 
formula 

= p[ x * - s t 7(x*)] (15) 

has the following property: either an index t* exists such that x 1 * G D*, or 
lim t ^ +00 f l best - f* = 0, where 

fLst = min /(x*) and f* = mmf(x). 

i=l,. ..,t xGD 

Let us check that g can be minimized through the projected subgradient 
method. First of all, g is convex with mhiagA™,! <?(«) = v(T, w) (see [7]) and 
we can easily show that u a G d a g(a), with u a bounded. In fact, for each 



point a G A m _ x the vector u a satisfies (12): 



g(a) — g(a) = max 2, a j^{Bj) — (u a , a) 
jeM 



jeM 

= (u a , a) — (u a , a) 
= (u a , a — a). 
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We now adapt the general updating rule (15) to our situation. For any 
diminishing step size rule {s t }fj^ and any vector a 1 G A m _i of coefficients, 
the update rule becomes 

a t+1 = p [ a * - stu 1 } = (a 1 - s t u l + Al)+ (16) 

where A e R is the normalizing constant such that 

rn 

£)(a*-s t «? + A)+ = l. 

i=l 

o 

Suppose now that a* GAm-i an d that the step size s t is sufficiently small to 
guarantee 

a\ - s t u\ + A > for each i G M. (17) 
Here A has to be chosen so that 

m 

^ - StUt i + A ) = X > 

i=l 

i.e., 

m m 

- s t ^2,u\ + m\ = 1, 

i=l i=l 

hence 

A = s t u\ 

where the average of the subgradient vector components. In 

what follows we will make sure to choose a diminshing step size rule small 
enough so that (17) is verified, or, equivalently, 

a\ - s t {u\ — u*) > for all t G N and for alii = 1, . . . , m. (18) 

We are now able to state the first convergence result. 

Proposition 4.5. Suppose (CD) and (MAC) holds. Let AJ n _ 1 be the bounded 
set of minimum points for g and let {s t }fj£ be a diminishing step size rule. 
Then, there exists another diminishing step size rule s t < s t which satisfies 
(18). Consequently, given a G A m _i and the recursive sequence 

it* = PVV{a t ) 

*+! = ot - s' t {u l - «*) (19) 



a 



either a G A* n _ x for some t* G N ; or 

lim a 1 = a* G AJ rt _ 1 and lim g(a t ) = g(a*) = v(T, w). 

t— ¥+OC " t— > + OD 



12 



Proof. First of all, notice that constraint (18) involves only those indexes 
i G {1, . . . , m} for which u\ > u l . Let be X t the set of those indexes in the 
step t. For each of them we would get 



OL; 



ma. 



uj-tt* (m - l)u\ - J2j^i ut j 
Now, fix an arbitrary integer K eN and define 

K-l 



ma. 



Ti 



mm 



K iei t I (m-lK-E^Wj 



Hence, let us define 

s' t = mm{s t ,T t }. 

Thus, s' t satisfies (18) and (13), since \im t ^ +00 s' t = lim 4 _> +00 s t = 0. 
To show ( fT4| ) , let us suppose Ylt=o s 't < +°°- This implies for some i* G 
{1, . . . , m} and some sequence {t r } C N 



lim 



maZ 



0. 



Since (m — l)uf — u / r > 0' taking a further subsequence {t p } C {t r } 

we have 



a-* — )■ 5j* = 0, so a* 3 ' — >■ a, = 1, and 

— >■ Uj for all j G {1, . . . , m}, 
with (m — l)£tj. > Mj. 

By continuity, w = (tti, . . . ,u m ) lies on the upper surface of V, so 

5>,->i. 

Moreover, w is supported by the hyperplane X]j=i = ^- 
By (20) and (21) we have that Ui* > — > 0, since 



m 



l)Wj» > Uj = Mj — Ui* > 1 — Mj. , so 



(20) 



(21) 



u^ > — > 0. 
m 
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Now, the coexistence of Ui* > and a>i* = clashes with the hypoth- 
esis (MAC). In fact, (u 1 ,...,u m ) G argmax^p J2jeM ®j x v m n j'\, = 
YljjLi* ^jUj = k and there is no (xi, . . . , x m ) G "P for which ^- gM > ^ 
Since > — > 0, there exists A? such that = > -j-. By 

(MAC) we can derive a partition from Ai* of (m — 1) subsets {-Sj}^*, with 
Uj^Bj = A { * and £?, n 5; = if j ^ /, such that fij(Bj) > e > for all 

If we consider the partition A defined as Ai* = 0, Aj = Aj U Bj, we get 
Yl & M A) = J] + ^-(-Bj)) = A; + (m — l)e > fc, 

which is a contradiction. 



The final statement is a direct consequence of Propositions 3.3, 4.4 and 
of the fact that if s' t — > then the sequence {a*} must converge to some a* G 
Aj n _ 1 . Moreover, by the dominated convergence theorem, g(a l ) — > g(a*) = 



v(T,w), last equality being again a consequence of Proposition 3.3, □ 



To prove the convergence of the PVVs and of the lower bound, we assume 
relative disagreement (RD). 

Proposition 4.6. Suppose (CD), (MAC) and (RD) hold. Then, for any 
a G A m _x and the recursive sequence (19), one of the following two condi- 
tions hold: 

• either u 1 * = (v(T,w), . . . ,v(T,w)) for some t* G N and v_{u t ") = 
v(T,w), 

• or lim t _>. + oo u* = (v(T, w), v(T, w)) and yju 1 ) —¥ v(T,w). 

Proof. By (RD) we have that for any point on Pareto border of V there 
exists one and only one hyperplane touching V (see [2]). By the conclusions 



of Proposition 4.5 either a 1 * G A^-i or lim t ^ +OC) a 1 = a* G A 



m— 1 ' 



In the first case, there exists only one PVV u l * corresponding to a 1 . Since 
the hyperplane with coefficients vector a 1 * touches the partition range V in 
the point corresponding to the maxmin allocation, then u l * must coincide 
with 

u* = (u(r, «>),..., v(I>)). (22) 

Also v(u*) = v(T,w). In fact, all the coordinates of u* are equal and, there- 
fore, maximal. Without loss of generality we choose u\ as maximal, and 

v(u*) = Ul ui _ - = Ul = v(T,w) 
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the last equality holding by ( 22 ) . 

In the second occurrence, suppose on the contrary that a 1 * — > a* while 
u f u*. Since the sequence {u 1 } is in a compact set, there must be a 
converging subsequence u l —> u ^ u* . The vector u is a second PVV asso- 
ciated to a*, but this is ruled out by (RD). Thus, \\m t ^ +00 u t = u*, and, by 
continuity, \im t _ ¥+oc v(u t ) = y(u*) = u*. □ 



4.1 The algorithm 

We now present two versions of an algorithm for the maxmin division prob- 
lem. The common initializing elements for both versions are listed in Table [Tj 
The first version computes upper and lower bounds for v(T,w) and updates 
the coefficient vector a through the subgradient rule p6| ). Both bounds are 
updated by means of a simple comparison with the old ones. The generic 
step is described in Table [2] A simpler but slower version, described in Table 
[3j computes the approximating optimal partition as well as the value. The 
finiteness of both algorithms is guaranteed by Propositions |4.4[ |4.5| and 4J3 



Particular care is needed in choosing the step sequence {s t }^. A se- 
quence converging too fast to may lead to an increase in the number of 
steps needed, since the step may soon become too small to reach the opti- 
mum. Similarly a sequence converging too slowly may result in values of a 
(and of the corresponding PVV's) jumping from one extreme to the other 
of the unit simplex (and of the partition range for the PVV's). This, again, 
will slow the convergence process. 



Table 1: Description and initialization algorithms elements. 



Elements 


Description 


Initialization 


a 


supporting hyperplane coefficients vector 


a° = l/m, j G M 


u 


PVV vector associated to a 


a 
U = U 


ub 


upper bound 


g(oP) 


lb 


lower bound 


v(u°) 


t 

s t 


diminishing step size rule 


s' = min{s , r } 



4.2 A five players example 

Let us consider the coalitional game defined in Q, with five players and 
players' preferences listed as probability distributions on C = [0, 1] in Table 
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Table 2: Generic step t for algorithm returning v(T,w). 


Step t 


Computation 


update a 


a* = a' -1 — ^(V -1 — u 1 ^ 1 ) 


update u 


yt = PVV^a 1 ) 


update ub 


if g(a t ) < ub, then ub = g(a t ) 


update lb 


if t>(w*) > lb, then lb = v_{u l ) 


stop condition 


if ub — lb < e, then STOP. Else, repeat the step. 


Table 3: Generic step t for algorithm returning u*. 


Step t 


Computation 


update a 




update u 


u * = pyy( a *) 


stop condition 


if maxjgjvf w* — minjgMU* < e, then STOP. Else, repeat 




the step. 



In Figure |2j we represent the initial densities (a) and then the maxmin 
partition for the fully competitive context (b), where Y = {{i}^^} and 
w({i}) = 1, for all i e N. 




0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



(a) Densities of preferences. (b) Maxmin partition. 

Figure 2: Densities of players preferences and maxmin partition in competitive context 
(red for player 1, green for player 2, blue for player 3, brown for player 4, purple for player 
5). 

For any S C N, we run our algorithm enforcing the two weight systems 
w car d and w pre , with a tolerance of 10~ 3 and we compute the corresponding 
game values (Table [5]). Consequently, in Table |6j we compute the Shapley 
value for each game. 
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Table 4: Players preferences. 



Player {i} 


Hi 


1 


Beta(2,5) 


2 


Beta(3,8) 


3 


Beta(7, 2) 


4 


Beta(10,10) 


5 


Uniform[0, 1] 



Table 5: Comparison between the coalitional game values. 



Coalition S 


77 (S, w card ) 


r)(S,W pre ) 


tt\ 7 — 1 2 S 4 5 


404 


0, 404 


{1,2} 


0,822 


0,842 


{1,3} 


0,835 


0,836 


{1,4} 


0,844 


0,861 


{1,5} 


0,819 


0,827 


{2,3} 


0,820 


0,820 


{2,4} 


0,826 


0,826 


{2,5} 


0,828 


0,833 


{3,4} 


0,808 


0,808 


{3,5} 


0,926 


1,040 


{4,5} 


0,886 


1,004 


{1,2,3} 


1,262 


1,280 


{1,2,4} 


1,273 


1,302 


{1,2,5} 


1,256 


1,265 


{1,3,4} 


1,275 


1,289 


{1,3,5} 


1,392 


1,465 


{1,4,5} 


1,366 


1,427 


{2,3,4} 


1,242 


1,241 


{2,3,5} 


1,389 


1,474 


{2,4,5} 


1,349 


1,414 


{3,4,5} 


1,403 


1,625 


{1,2,3,4} 


1,706 


1,727 


{1,2,3,5} 


1,877 


1,903 


{1,2,4,5} 


1,841 


1,862 


{1,3,4,5} 


1,968 


2,044 


{2,3,4,5} 


1,940 


2,032 


{1,2,3,4,5} 


2,477 


2,477 



The two games share the same ranking for the Shapley values 
P1.5 y P1.3 y P1.4 >- Pl.l y P1.2 
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Table 6: Comparison between the Shapley values. 



Playerji} 


S.V. in 7](S,w card ) 


S.V. in rj(S,w pre ) 


1 


0,465 


0,436 


2 


0,451 


0,425 


3 


0,507 


0,519 


4 


0,491 


0,502 


5 


0,563 


0,594 



which therefore seems to be robust enough to the choice of system weights. 

Also, the weight system w pre amplifies the difference in the Shapley values 
obtained with w car d, yielding a higher variance for the values' distributions. 

5 Concluding remarks 

In the previous section we described a couple of algorithms that return 
maxmin values and partitions in both competitive and cooperative settings. 
It is important to note that we could think of the same procedures as interac- 
tively implemented between (coalitions of) players and an impartial referee. 
At first the referee proposes a division of the cake based on the maxsum 
division of the cake with equal weights for all players. The players now re- 
port their utilities and the referee corrects the inequalities in the division by 
proposing a new maxsum division with modified weights: Players who were 
better off will be given a smaller weight and those who were worst off will see 
their weight increase. Of course, one cannot hope to achieve the same degree 
of precision, since the algorithm performs that step dozens of times, but the 
bounds described in Section 3 give a precise idea on how far the proposed 
division is from the desired one. 

Many issues remain open. We hint at two of them. 

• In the numerical example it would be interesting to link the Shap- 
ley value rankings to the original system of preferences. What makes 
Players 5 and 3 the most powerful players in the cooperative division 
process? Apparently the two utility functions have different features: 
Player 5's distribution is uniform over the unit interval and his density 
is maximal only at the very ends of the interval. On the other hand, 
Player 3's preferences are concentrated at the second half of the inter- 
val - where he has no competitors, except player 5 (who, however, has 
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a smaller density). No simple explanation could be provided so far. 

• Beyond the convergence of the algorithms, which end in a finite number 
of steps, returning the approximate solution up to a specified degree 
of precision, it would be interesting to investigate about the computa- 
tional efficiency of the same algorithms 
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